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Abstract 

We have studied the analytical Markovian and non-Markovian dynamics of 
quantum correlations, such as entanglement, quantum discord and Bell non- 
localities for three noisy qubits. Quantum correlation as measured by quan- 
tum discord is found to be immune to death contrary to entanglement and 
Bell nonlocality for initial GHZ- or W-type mixed states. 
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1. Introduction 

Entanglement, quantum discord and Bell nonlocalities and their dynamics 
under the influence of external noise have been an active area of research, 
recently 0, i, i, S, H 0, S, S 0, [nl, H, El . Although entanglement was 



considered the most important form of quantum correlations for carrying 
out quantum computations [13] or other quantum tasks, such as quantum 



teleportation [15| and quantum cryptography [16[, it was shown by Datta et 



al. that states that have no entanglement but non-zero quantum discord still 



can be used to perform useful quantum tasks [171 . 

Entanglement, quantum discord and Bell nonlocality of a given state are 
all different aspects of quantum correlations in that state. Entanglement 
refers to separability of the states while quantum discord is the difference 
between the total and the classical correlations for a given state as measured 
by von-Neumann entropy. Quantum discord is an independent measure of 
nonclassical correlations which might include entanglement, but the rela- 
tion between entanglement and quantum discord is complicated for a general 
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states 



18l . Il9l . I20I |21| . Bell nonlocality is involved with the violations of Bell 



inequalities which identify the entangled mixed states whose correlations can 



be reproduced by a classical local model |22 



All quantum systems have an environment that entangles with the system 
and drives it to an incoherent state. It was shown by Yu and Eberly that while 
coherence of a single qubit goes to zero exponentially, the entanglement be- 
tween two such qubits can cease to exist in a finite time, a phenomenon called 
entanglement sudden death (ESD) [2^ 



23 , a number of studies 



After Ref. 

on the dynamics of other types of quantum correlations, such as quantum 
discord and Bell nonlocalities under Markovian as well as non-Markovian 
memory types were carried out to investigate whether sudden death phe- 
nomena exist for those correlations ^ [a, [s], [q], 0, 11, 12, 13 1. Wang et 
al. have considered the effects of a non-Markovian dissipative environment 
on the dynamics of a two qubit system for Bell-like or extended Werner-like 
initial states and found that quantum discord is more robust compared to 
the entanglement as measured by concurrence 10|. Werlang et al. stud- 



ied the dynamics of quantum discord and entanglement for the Markovian 
environment with depolarizing, dephasing and generalized amplitude damp- 
ing channels and found that whenever there is ESD, the discord decreases 



only exponentially [ll| . A number of work dealt with the comparision of the 



behavior of entariglement and Bell nonlocality for two qubits in various en- 



vironments 



0,0: 



12, 22 



Werner has demonstrated that mixed-entangled 
states might not violate Bell inequalities [22[. Bellomo et al. have found 
that Bell inequality might not be violated for a state with high values of 
entanglement for a two qubit system subject to amplitude damping 0]. 

Although the relation between the memory properties (Markovian and 
non-Markovian) and type of correlation and the resulting dynamics is inves- 
tigated for a bipartite systems by a large number of groups, the works on 
multipartite systems are scarce [l|, [i], li], 0, 0] • Among them, Ann and Jaeger 
examined the finite time destruction of multipartite Bell inequality violation 
in three-qubit systems initially prepared in GHZ and W pure states under the 
multi-local asymptotic dephasing noise in Markovian approximation [l|, lij. 
Moreover, Qui et al. showed that the tripartite Bell inequality violations 
can be fully destroyed in a finite time for a general pure state in a nonlo- 
cal antiferromagnetic environment under an external magnetic field by using 
spin wave approximation [3}. Also, Yang et al. demonstrated the three qubit 
Bell nonlocality sudden death in Tavis-Cummings model with cavityloss for 
W-like initial state by solving the Schrodinger equation exactly Liu, 
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Shao and Zou considered the dynamics of Bell-nonlocality for a three-qubit 
Heisenberg XY chain in the presence of an external magnetic field for pure 
GHZ-state and found that the multipartite Bell-inequality violations can be 
destroyed in a finite time because of intrinsic decoherence [sj. However, none 
of them compares the dynamics of quantum correlations. 

In the present study, we have analyzed the dynamics of Bell nonlocalities 
as measured by Mermin-Ardehali-Belinksii-Klyshko (MABK) and Svetlichny 



inequalities |24j. |25| . |26| . l27l |. quantum discord and entanglement as measured 
by concurrence for bipartitions and tripartite negativity for tripartite states 
for a system of three qubits which have energy levels that are stochastic with 
Ornstein-Uhlenbeck type correlations. The analytic expressions for quantum 
discord. Bell nonlocalities as well as concurrence and tripartite negativity are 
derived for W- and GHZ-type initial states by exploiting a procedure based 
on the knowledge of single-qubit dynamics j28|. We consider both Markovian 
and non-Markovian time evolution and compare and contrast the character 
of time dependence of quantum correlations. 

The organization of this paper is as follows. In Sec. 121 we introduce the 



model and its solution using the procedure analyzed in Appendix Appendix A 
In Sec. [31 we show explicit analytic calculations of quantum discord, con- 
currence, tripartite negativity and Bell nonlocalities for three-qubit system 
initially prepared in GHZ- and W-type states and explore the effects of non- 
Markovianity and mixedness on quantum correlations. We conclude with a 
summary of important findings in Sec. |H 

2. The Model and its Solution 

The model we consider in this paper is three uncoupled qubits interacting 
with their noisy environments independently which cause time-dependent 
fluctuations in their energy levels. The Hamiltonian for this model is given 



by |29| (we set h = 1): 



where dz is the usual Pauli spin operator in z-direction, flA,B,c{t) are the 
independent fluctuations of the transition frequencies obeying non-Markovian 
approximation with mean value properties 

M{a(t)} = 0, (2) 
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M{ni{t)Qi{s)} = a{t-s) 

= ^e-'^\'-^\,^ = A,B,C, (3) 

where M{...} stands for the statistical mean over the noises flA{t), ^B{t) 
and flcit). Here Ti{i = A,B,C) are the damping rates due to the couphng 
to the environments, 7 is the noise bandwidth which determines the finite 
correlation time of the environmental noise (r^ = 7~^) and a{t — s) is the 
reservoir correlation function. For simplicity, we will take the noise properties 
to be the same for A, B and C (e.g., T a = ^ b = Tc = T)- And note that 
in the limit 7 —)■ 00 (tc —)■ 0), Ornstein-Uhlenbeck noise reduces to the well- 
known Markovian case 

a{t- s) = V5{t- s). (4) 

For the total system described by the Hamiltonian ([1]), the stochastic Schrodinger 
equation is given by 

^J^\^{t))=H{t)\^>{t)), (5) 

with solution 

|^(t)) = ?7(t,fiA,^^B,^^c)|^(0)), (6) 
where the stochastic propagator f/(t, VLa, ^b, ^c) is given by 

u{t,nA,nB,nc) = e-t/o(f^AW^^+f^Bw^f+now^?)d._ (7) 

The reduced density matrix for spins A, B and C is then obtained from the 
statistical mean 

m = Mmt))m)\}- (») 

Since we assume independent reservoirs and the qubits interact only with 
their own independent reservoir, we can use the procedure developed by Bel- 
loma et al. to obtain the time dependent density matrix of three qubit system 



from the knowledge of the dynamics of single qubit in an environment [28 



Using the stochastic Schrodinger equation ([5]), the master equation for the 
single qubit reduced density matrix can be derived as [3l|, l32|, l33[ : 



|pW = -^(pW-^.pW^.), (9) 
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where 

r 



G{t)= I a{t-s)ds = ^{l-e-^"). (10) 







Here G{t) is a time-dependent coefficient which includes the memory infor- 
mation of the environmental noise. The analytic solution of Eq. ([9]) and the 
time-dependent density matrix for tripartite system obtained according to 



procedure of Ref. [28|] is outlined in Appendix Appendix A 

Since we assume the noise properties to be the same for all qubits, the 
parameters in Eq. A flA.Sp can be obtained from the solution of Eq. with 
the help of Eqs. A dAlTl) and A dAlSl) as: 

uf = uf = uf' = l, 

vf = vf=v^ = 0, 

zf = .f = .f = e-/W, (11) 

where 



fit) = fG{s)ds 
Jo 



^ ^t+-{e-^'-l)]. (12) 
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Using Eqs. A (1A.9I) . A llA.lOl) and (ITTl) . the matrix elements of the reduced 
density matrix for three-qubit can be determined: 

Piiit) = Pu{0), {z = l,2,...,8}, 

Pij{t) = pij{0)e~^^*\ {ij = 12,13,15,24,26,34,37,48,56,57,68,78}, 

p.^.(t) = pij{0)e-^f^^\ {ij = 14,16,17,23,25,28,35,38,46,47,58,67}, 

Pijit) = p.,(0)e-3^W, {zj = 18,27,36,45}. (13) 

3. Correlation Measures 

In this section, we briefly review the correlation measures considered in 
the present study. 

Quantum Entanglement: For a pair of qubits, the concurrence as a 



measure of entanglement is well-defined. It is introduced by Wootters [34 
and defined as 



G = max{0, VAi - VA2 - VA3 - VA4}, (14) 



5 



where the quantities Aj(i = 1, 2, 3, 4) are the eigenvalues of the matrix 



ABK^y 



(15) 



in descending order and p^^ is the conjugate of pab- For the density matrix 



with X-form 29 



Pab 



a w 

6 z 

2 c 

w Q Q d 



(16) 



the concurrence function f ll4p has a simple analytic form: 



C = 2 max{0, z 



ad, w 



be}. 



(17) 



On the other hand, for the three-qubit case concurrence is not well- 
defined. Its calculation is based on numerical optimization procedure which 
does not guarantee exact results |35|, l36| . Instead, for mixed states of three 
qubit systems, the tripartite negativity N was introduced by Sabin and 
Garcia-Alcaine l37| as 



A-BC^B-AC^^'C-AB 



Nr 



,1/3 



where Nj^jk with / = A, B,C and JK = BC, AC, AB is the negativity 
with respect to the subsystem / and defined as [38 



I-JK 



(19) 



where ai{ff^) are the negative eigenvalues of the partial transpose p*^ of the 



system density matrix with respect to the subsystem / [39 



One should note that for the mixed states, N is not able to quantify 
multipartite entanglement fully, but its positivity ensures that the state under 
consideration is not separable [40|. N is used to study the time evolution 
of entanglement for a number of tripartite systems 40|, |4l|, |42|, |43|. For the 



type of GHZ and W states and the dynamics considered in the present work 
is the same as Nj^jk 



43 



Quantum Discord: For two-qubit systems, quantum discord as a mea- 
sure of quantum correlation was introduced by Ollivier and Zurek fisl]. It 
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is defined as the difference between two expressions of mutual information: 
total and classical correlations, namely, 

D = I{pab)~J{pab)- (20) 
Here I{pab) is the total correlation between two subsystems defined as 

Hpab) = S{pa) + S{pb) - S{pab), (21) 

where S{p) = —Tr{p\og2p) is the von Neumann entropy and pa{pb) is the 
reduced density matrix of pab obtained by tracing out B{A) [39]. The other 
quantity J{pab) is the classical correlation between A and B as the maximum 
information one can get from A by measuring B. It is defined as 

J{Pab) = max{S{pA) - ^^^^(p^iJ}, (22) 
{n.} V 

where pa\ = r!^^'^^\k^^^^^\ and flj is a set of projectors defined as flj = 

J(g)|i) {i\ (i = 1,2) where |1) = cos6' |+^) + e*'^sin6' j-^) and |2) = sin6' |+^)- 
e*'^cos6' |— 2) are the orthogonal states [9]. The projector operators are used 
to measure subsystem B, corresponding to the outcome i with probability 

Pi = TrABiJiiPABfii)- 

In general, it is hard to calculate the analytic expression for the quantum 
discord. However, if the reduced density matrix of two qubits has X-form as 
Eq. (fT6|l with b = c, the quantum discord has a simple analytic form 



D = mm{Di,D2}, (23) 



where 



Di = S{pa) - S{pab) - a'^og2 [ — -61og; 



a + b J \a + b 



d^og2 [ r—l) -b\og2 



b + dj \d+b 
D2 = S {p a) - S {p ab)- 1(1 + ^^)log2(i- {I + k] 



2' ' \2 
where = {a — d)"^ + 4(|z| + \w\) 



^l-«:)logJi(l-«:)), (24) 
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Bell-Nonlocality: For three-qubit case, violation of Bell nonlocality can 
be used as a measure of quantum correlations. From this point, we choose 
two kinds of nonlocality measures. The first one is introduced by Mermin- 
Ardehali-Belinksii-Klyshako (MARK) [l|, S, H [isl |26|- It is easily com- 



putable for three-qubit case and violated whenever 
1 where MABK operator is given by 
1 



B 



TT{Bp{t)) 



> 



B = - [ MaMbM'c + MaM'bMc + M'^MbMc - M'^M'^Ad'c 



(25) 



The second inequality is put forward by Svetlichny [H [i], 27| and denotes 



Tr{Sp{t)) 



> 4 where the 



genuine tripartite Bell nonlocality if 
Svetlichny operator is 

S = MaMbMc + MaMbM'c + MaM'bMc + M'^MbMc 

- M'j^M'bM'c - M'aM'bMc - M'j^MbM'c - MaM'bM^^ (26) 

where Mk is the measurement operator for the K^^ qubit, the primed and 
unprimed terms correspond to different measurement directions for the mea- 
suring party. The measurement operator for each successive subsystem is 
obtained from the preceding one by a rotation: 



Mk 
M'k 



cos 6 k 
sin Ok 



- sin 6k 
cos Ok 



Ma 
M'a 



(27) 



where 6k {K = B, C) are the rotation angles [H, [ij. 

To demonstrate the finite time loss of Bell nonlocality in both Markovian 
and non-Markovian environments, we shall assume two different nonlocality 
regimes satisfying: 



>4, 



<4, 



where the regime (i) shows that the state p(t) is genuinely tripartite Bell 











> 1, 











ii) 




< 1, 













nonlocal. Also we choose the rotation angles in which 
have maximum values at t = 0. 



B 



and 



S 
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In the following, we will consider the dynamics of tripartite GHZ- and 
W-type states. As is well known, GHZ-state cannot be transformed into 
W-state by local operations and classical communication 4J]. This means 
they bear incompatible multipartite correlations. From this point of view, 
one may expect that the two initial states might have some differences in 
their correlation dynamics due to effect of external noise. 

3.1. GHZ-type initial states 

In this section, we explore the effects of noise on the following three qubit 
GHZ-type initial states, 

1 — r 



m 



-h + r\GHZ) {GHZ\, 



(28) 



where \GHZ) = (|111) + |000)) is the GHZ-state, r is the purity which 

ranges from to 1 and Is is the identity matrix of dimension 8. For this 
initial state one can define the measurement operators as Ma = o"j, and 



1 , then 



Ma 


= ay 


® /2 « 


/2, 




M'a 


= dr^ 




h, 




Mb 


= 12 


® [cos( 


dB)cry - 


■ Sm{9B)dr:j] ® 12 


M'b 


= h 


® [sin( 


0B)o-y + 


cos(6'b)o-^] (8) /2 


Mc 


= 12 


® /2 ® 


[cos{9c 


)cry - sm{9c)cr^] 


M'c 


= h 


® /2 ® 




&y + cos(6'c)o-^.] 



(29) 

The expectation value of B and S operators for the GHZ-state under the 
classical noise can be easily obtained as: 



B 



S 



2re-'f^'^\cos{9Bc)\, 
4re-^^(*) |cos(0Bc) - sin(^Bc) 



(30) 



where 9bc = 9b + 9c- Also, it should be noted that this state does not 
have any bipartite correlations, however, it has non-zero tripartite negativity 
which is equal to 



1 



N = max{0, — (1 - r - 4re"^-^(*))}. 



(31) 
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Figure 1: The dynamics of iV (solid plots), 



B 



- 1 (dashed plots) and 



-4 (dotted 



plots) versus Yt with r — 0.98 for GHZ-type initial state. Here the thick plots correspond 
to non-Markovian regime with 7/r = 0.1 and the thin plots to Markovian regime with 
7/r = 10. 




Figure 2: The dynamics of 



B 



1 (Fig. (a)), 



4 (Fig. (b)) and N (Fig. (c)) 



versus Ft and Markovianity, 7/r, with r ~ 0.98 for GHZ-type initial state. 



In Fig. 1 and 2, we have displayed the effects of non-Markovianity on 
the dynamics of tripartite entanglement and Bell nonlocalities for GHZ-type 
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initial state with r = 0.98. Both figures indicate that tripartite entangle- 
ment and Bell nonlocality violations suffer sudden death. It can be seen 
from Fig. 1 and 2 that the lifetime of tripartite entanglement as measured 
by the tripartite negativity is significantly longer than the lifetime of Bell 
nonlocalities in Markovian as well as non-Markovian coupling regimes. The 
non-Markovianity only delays the death, not prevent it, as can be seen from 
Fig. 2. 




Figure 3: The dynamics of (B 

and N (Fig. (e) and (f)) ver 
correspond to Markovian regime with 7/r 
regime with 7/r — 0.1. 



1 (Fig. (a) and (b)), 



4 (Fig. (c) and (d)) 



sus Tt and r for GHZ-type initial state. Fig. (a), (c) and (e) 
10 and Fig. (b), (d) and (f) to non-Markovian 



The effect of purity is considered in Fig. 3, where we have plotted the 
tripartite negativity and Bell nonlocalities versus dimensionless time and the 
purity, r, for GHZ-type initial state in Markovian (7/r = 10) and non- 
Markovian (7/r = 0.1) coupling regimes. It can be seen that the purity is 
the most important parameter for the existence of Bell nonlocalities and the 
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entanglement. 











and 











suffer death at all time points for r < 0.5 



and r < 0.4, respectively, while the non-zero values of tripartite entanglement 
exists for a much wider range of purity values, 0.2 < r < 1. The range of r for 
nonzero and Bell-inequality violations is independent of the Markovianity 
of the dynamics. Another important point is that for pure state (r = 1), the 
Bell nonlocalities cease to exist in a finite time while entanglement dies only 
asymptotically. 

3.2. W-type initial states 

The W-type initial state can be expressed as: 



m 



1 — r 



is + r\W) {W\ 



(32) 



where \W) = -^(|100) + |010) + |001)) is the W-state and r is the purity. For 
this type of initial states, the measurement operators for Bell nonlocalities 
are given by j2| 



Ma 




5 /2 (g) h, 


M'a 


= cr^ C 


/2 ® h, 


Mb 


= /2$ 


5 [cos(6'b)(T^ - sin{6B)cra:] (8) h, 


M'b 


= /2$ 


5 [sin(6's)(T^ + cos(6'B)a-^] (g) /a. 


Mc 


= /2$ 


5 h (g) [cos(6'c)(3-2 - sin(6'c)<3-x], 


M'c 


= /2$ 


5 I2 ® [sm{6c)crz + cos{6c)ax]. 



(33) 



Then the time- dependent expectation values of the operators B and S can 
be calculated as: 



B 



S 



^l + 2e-2/W) |sin( 



^BC) 



;i + 2e-2/W) |(cos(0Bc) + sm{9Bc))\ ■ 



(34) 



Unlike GHZ state, W state has a high degree of bipartite correlation. 
Thus, W-type states allow some comparison between tri- and bi-partite cor- 
relations. The analytic forms of the tri- and bi-partite entanglements and 
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quantum discord for this state can be easily calculated as: 



N = — max{0, -3 + 3r + 8y2re~^-^(*)}, 

C = ^max{0,4re"2-^(*^ - ^3(1 - r)(3 + r)}, 
6 

D = mm{Di,D2}, 



(35) 



where 
Di = 



i=l 1=1 ^ 



a + J \ a + 



dlog2 



d 



b + d 



61og2 



d + b 



E l°g2(Af ) + Xt^ hg.iXn -1(1 + ^) lo& (^(1 + n)) 

1=1 i=l ^ ^ 



i(l - k) log2 (^(1 - ^ 



3+r \AB 
12 ' ^3 12 



Tl(3 + r - 

2/(i) 



where A^^ = 3^,A^ = ^.A^^^ = if^.A^^ = 
4re-2/W),Af^ = ^(3 + r + 4re-2/W), a = if:, 6 
and K = |Vl + 4e-4/W. 

The Markovian and non-Markovian dynamics of Bell nonlocalities, bl- 
and tri-partite entanglement and quantum discord for W-type initial state 
are displayed in Figs. 4(a) and (b), respectively for r = 0.98. One obser- 
vation from these figures is the fact that tripartite entanglement as mea- 
sured by tripartite negativity has a longer lifetime compared to the bipartite 
entanglement as measured by concurrence, independent of Markovianity of 
the dynamics. One can also deduce from Figs. 4(a) and (b) that both Bell 



nonlocalities 



B 



and 



S 



are more fragile than all the other quan- 



tum correlations considered in the present study. Also, quantum discord is 
found to be immune to sudden death independent of the Markovianity of 
the dynamics. One should note that quantum discord has a discontinuity 
at rt = 0.4 for Markovian and Tt = 2.3 for non-Markovian dynamics. The 
reason for the discontinuity is in the definition of quantum discord which 
involves min{Di, D2} where Di and D2 are defined in the text. Ref. 45] has 
shown experimentally that the discontinuity is real. 
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Figure 4: (Color online) The dynamics of N (thick blue solid plots), C (thin black dotted 



plots), D (thin red solid plots) 



B 



1 (thick dashed red plots) and 



S 



4 (thick 



black solid plots) versus Tt with r = 0.98 for W-type initial state. Here Fig. (a) corresponds 
to Markovian regime with 7/r = 10 and Fig. (b) to non-Markovian regime with 7/r — 0.1. 




Figure 5: The dynamics of 



B 



1 (Fig. (a)) and 



7/r with r = 0.98 for W-type initial state. 



S 



4 (Fig. (b)) versus Ft and 



The effect of Markovianity on Bell nonlocalities and tripartite negativity, 
concurrence and quantum discord for W-type initial states are displayed as 
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Figure 6: The dynamics of N (Fig. (a)), C (Fig. (b)) and D (Fig. (c)) versus Ft and 7/r 
with r — 0.98 for W-type initial state. 



functions of Tt and 7/r in Fig. 5 and 6, respectively. The effect is found to be 
a prolongation of lifetime for the correlations that suffer sudden death (Bell 
nonlocalities, and C). 




Figure 7: The dynamics of 



B 



-1 (Fig. (a) and (b)) and 



S 



4 (Fig. (c) and (d)) 



versus Ft and r for W-type initial state. Fig. (a) and (c) correspond to Markovian regime 
with 7/r = 10 and Fig. (b) and (d) to non-Markovian regime with 7/r = 0.1. 
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We have analyzed the purity dependence of Bell-nonlocahties as measured 



by 



B 



and 



S 



for the Markovian and non-Markovian dynamics and 



present the resuhs in Fig. 7 (a)-(d). One can conclude two important find- 
ings from these figures: i) Both of nonlocalities are highly purity dependent; 



B 



suffers death at all times in the regions where r < 0.7, while the 



non-zero values of 



S 



is limited to a much narrower range of r values 



(0.9 < r < 1). These are independent of Markovianity of the dynamics, ii) As 
expected, non-Markovian dynamics offer a longer lived nonlocality compared 
to the Markovian case. 




Figure 8: The dynamics of N (Fig. (a) and (b)), C (Fig. (c) and (d)) and D (Fig. (e) 
and (f)) versus Ft and r for W-type initial state. Fig. (a), (c) and (e) correspond to 
Markovian regime with 7/r = 10 and Fig. (b), (d) and (f) to non-Markovian regime with 
7/r ^ 0.1. 
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Similar to the case of Bell-nonlocalities, we have considered the purity 
dependence of Markovian and non-Markovian dynamics of tripartite nega- 
tivity, bipartite concurrence and quantum discord and displayed the results 
in Fig. 8 (a)-(f). The most important finding from Fig. 8 (a)-(f) is the nonex- 
istence of death for the quantum discord independent of purity of the initial 
states and Markovianity of the dynamics, which indicates that although en- 
tanglement type quantum correlations die out in finite time, not all types of 
quantum correlations are lost because of the noisy environment. Ferraro et al. 
showed that quantum discord can not have sudden death in Markovian envi- 



ronments 



13|. Our results indicate that absence of sudden death for quan- 
tum discord is valid for Markovian as well as non-Markovian dynamics for the 
particular states considered in the present work. Another important observa- 
tion concerns the behavior of Bell-nonlocalities and entanglement measures 

for the pure initial states (r = 1). While there is a sudden death of 



' p 



and 



for r = 1, there is no ESD in or C . Also, nonzero bipartite 



S , 

entanglement is found for a narrower range of initial purity (0.5 < r < 1) 
compared to nonzero tripartite entanglement (0.2 < r < 1). 



4. Conclusions 

We have analyzed the dynamics of quantum correlations, such as quan- 
tum discord, entanglement and Bell nonlocalities for three qubits that have 
stochastic time-dependent level spacings. The considered noise have Ornstein- 
Uhlenbeck type correlation. The dynamics is considered for GHZ- and W- 
type initial states and the survival times of the quantum correlations are 
compared. The tripartite entanglement is found to be immune to sudden 
death for pure GHZ as well as W states, while Bell inequalities cease to be 
violated for all types of initial states considered in this study. For GHZ- 
type initial states, there is no bipartite entanglement or quantum discord, 
while Bell nonlocality and tripartite entanglement as measured by tripartite 
negativity is nonzero for purity greater than approximately 0.4 and 0.2, re- 
spectively. Bell inequahty is found to be not violated at all purities at much 
shorter times compared to the lifetime of entanglement. 

W-type initial states display a richer dynamics, as they contain both 
bi- and tri-partite entanglement as well as Bell nonlocalities and nonzero 
quantum discord, initially. Quantum discord is observed to be more robust 
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compared to entanglement, because it decays exponentially while the concur- 
rence of the same state suffers sudden death. Also, tripartite entanglement is 
found to survive longer compared to the bipartite one if the state is mixed; for 
pure W states, both bi- and tri-partite entanglements decrease exponentially 
with time. 

We have considered the effects of non-Markovian dynamics of both types 
of states and found that the sole effect of non-Markovianity is to prolong 
the lifetime of quantum correlations compared to the Markovian dynamics, 
except for quantum discord which decreases only exponentially. Moreover, it 
should be noted that under Ornstein-Uhlenbeck type noise, once the entan- 
glement and Bell nonlocalities die, rebirth or revival does not occur in their 



dynamics [29 
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Appendix A. 

We consider a system tliat includes three subsystems 5* = A,B,C, and 
each subsystem includes one qubit, S = A,B,C, interacting with its local 
reservoir Rs and there is no other interaction in the whole system. Initially, 
each qubit and its reservoir are independent, thus the evolution of the reduced 
density matrix for qubit S is given by 

p'{t) = Trns{uHt)p'iO) ® p^%0)U~'\t)} , (A.l) 

where the trace is taken over the reservoir degrees of freedom and U^{t) 
is the time evolution operator for the part 5*. Eq. flA.ip can be expressed in 
terms of the Kraus operator K^j^ (t) as 

p'{t) = j2^!,{t)p'm'j,{t). (A2) 

a/3 

Since there is only interaction between single qubit and its corresponding 
reservoir, the time evolution operator U'^it) of the complete system factorizes 
as 

(/^(t) = U^it) ® U^it) ® U^{t), (A.3) 
then the reduced density matrix (1A.2P for one qubit can be extended for three 



qubit s as 

oi/3i a202 0:3/33 

Given the basis {\ls) , I'^s)} for qubit S, inserting the identity operator / = 
\is) {is\, one can get the reduced density matrix elements of the single 



qubit in the given basis as 

J si's 



{^s\p'\^'s)=pU,it) = (A.5) 



Isl's 



where = {^s\ K^^f^sit) \^s) (/5l^f]/3,W 14)- Combining Eqs. (0, (jAl 



and (lA.Sp . we finally get the reduced density matrix elements of the three- 
qubit system 

Pn.i..^.3.i it) = E E E (^)^S (^)^':| it)p,,,>Mi'2,isi', m-Q) 

ii^i ^3^3 
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The procedure given above allows us to obtain the dynamics of three 
qubits, provided that the dynamics of one qubit is known, by a purely alge- 
braic way and independently from the initial conditions. The procedure for 
N-qudit system is given in Ref. j28i |. 

We consider the single-qubit density matrix in the form, 

where 

pfiW = «fPn(0) + ^^f(t)pf2(0), 

Puit) = pll{t) = z^pUO), (A.8) 

where , and are functions of time and determined by the model of 
chosen. 

By considering Eqs. (]A.5|) . f lA.6|) and (]A.8I) . in the standard basis |1) = 
|111),|2) = |110),|3) = |101),|4) = |100),|5) = 1011), |6) = |010) , |7) = 
|001) , |8) = |000), we can obtain the diagonal elements of the reduced density 
matrix of three-qubit system as 

pii(t) = ufufufpniO) + ufufv^ p22i0) + ufvfufpssiO) + ufvfv^puiO) 
+ vfufufp^m + vfufv^peeiO) + vfvfufp77{0) + vfvfv^pssiO), 

P22{t) = ufuf{l-uf)pr,{0) + utuf{l-v^)p22{0) + utvf{l-uf)p,s{0) 
+ ufvfil - vf)p,M + vfufil - )p55(0) + vtuf{l - vf)p,,{^) 

P33{t) = Mf(l-«f)nfpn(0) + nf(l-Mf)t;fp22(0) + nf(l-i;,>fp33(0) 
+ uf{l - vf)v^p,,{0) + vf{l - uf)ufp,,{0) + vf{l - «f )t;f P66(0) 
+ vf{l - vf)ufprr{0) + vf{l - vf)v^pss{0), 

PM{t) = Ut{l-Uf){l-Uf)pii{0)+Ut{l-Uf){l-V^)p22{0) 

+ ufil - vf){l - )p33(0) + ufil - vf){l - v^)pu{0) 
+ vf{l - uf){l - uf)p,,{0) + vf{l - uf){l - v^)p,eiO) 
+ vf{l - vf){l - uf)pn{0) + vf{l - vf){l - )P88(0), 

P55(t) = {l-ut)ufufpu{0) + {l~ut)ufv^p22{0) + {l-ut)vfufp,s{0) 

+ (1 - uf)vfv^p,,{0) + (1 - vf)ufufp,,{0) + (1 - vf)ufv^p,e{0) 
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vf)vfufprm + (1 - vf)vfv^pss{0), 
ut)uf{l - uf)pn(0) + (1 - ut)uf{l - v?)p22(0) 

Ut)vf{l - «f )P33(0) + (1 - Ut)vf{l - V^)pu{0) 

vf)uf{l - «f )P55(0) + (1 - vf)uf{l - )P66(0) 
vf)vf{l - uf)p,r{0) + (1 - vf)vf{l - v^)pss{0), 
ut){l - «f Xpn(O) + (1 - ut){l - uf)vfp,M 
nf){l - vf)u'^p,m + (1 - uf){l - vf)v?p,M 
vt){l )«f p55(0) + (1 - vt){l - uf)vfp,e{0) 
vf){l vf)ufprm + (1 - - vf)vfpss{0), 

ut){l - uf){l - uf)pu{0) + (1 - ut){l - uf){l - vf)p,M 

Uf){l - Vf){l - Uf)pss{0) + (1 - Ut){l - Vf){l - )P44(0) 
Vf){l - Uf){l - «f )P55(0) + (1 - Vf){l - Uf){l - V^)p,,{0) 
Vf){l - V^){1 - Uf)p,,{0) + (1 - Vf){l - V^){1 - )P88(0), 

(A.9) 

and the off-diagonal elements as 

pu{t) = utufz^pniO) + ufvfz^puiO) + vfufz^pseiO) + vfvf p^s{0), 
Pisit) = u^z^ufp^-iiO) + ufzfv^p^M + vtzfu'^p^,{Q) + vtzfv^p^^{Q), 
pu{t) = utz^z^pu{0) + vfzfz^p,s{0), 

pi,it) = zfufufpi,iO) + zfufv^p26iO) + zfvfufp^riO) + zfv^vf p^siO), 
p^eit) = zfufz^p^e(0) + zfvfz^p,s{0), 

p,j(t) = Zfz^ufp,7(0)+zfzfv^p28(0), 

pisit) = zfzfzfprsiO), 

P2s(t) = ufzfzf^*p2m+vfzfz^*p,m, 

P24{t) = Ufzfil - )pi3(0) + Utz^il - V^)p24{0) + vfzf{l - uf)p,,{0) 

+ vfzf{l-v^)p,s{0), 

P2,{t) = Zfufz^*p2m+4vfz^'*p,j{0), 

P2e{t) = zfuf{l-uf)p,,{0) + zfuf{l-v^)p2e{0) + zfv^{l-uf)ps^{0) 
+ zfvf{l-v^)p,s{0), 

P27{t) = Ztzfzf*p2m. 

P2s{t) = zfzf{l-uf)p,r{0) + zfzf{l-v^)p2s{0), 

Pu{t) - U,^(l-Uf)^fpi2(0)+U,^(l-O^fp34(0)+^,^(1-O^fp56(0) 



Pm{t) 



P77{t) 



+ 

+ 
+ 
+ 

+ 
+ 



+ 
+ 
+ 
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P35(^) 

Psrit) 

pAhif) 
Pi&{t) 
P47{t) 
Pm{t) 

Pb&{t) 

P57{t) 

P58it) 
P67{t) 

Pesit) 

P78{t) 



+ 



+ 



+ 



+ 



+ 



+ vf{l-vr)zfprs{0), 



AB*„,C 



<P35(0) + ^f^r^;^P46(0), 

^^fp36(0), 



B\„,C' 



+ zf{l-vf)vf'p4s{Q), 



c . 



z^p,,{0) + zt{l-v^)z^pss{0), 

P45(0), 



(1 - «f )p35(o) + z^z^*{i - vi:')p,,io) 



c* 



p,,{0) + zf{l 



V, z 



c* 



P47(0), 



zfil - <)(1 - u^)p^m + zfil - <)(1 - t;np26(0) 
zfil - vf)il - uf)psm + zfil - vf)il - )P48(0), 



— u 



puiO) + (1 - O^f P34(0) + (1 - vf)ufz^'p,,iQ) 



vf)vfz^ P7si0), 



Vf)zfv^p,8i0), 



u 



^)zfzfp,M + {l-vf)zfz^p,m, 



B^C . 



f)zfzf*p,,iQ) + (1 - vf)zfzf*p,,iQ), 



1 - ut)z^v^p,,iQ) + (1 - ut)z^v^p2M + (1 - vt)z^u^p,,iQ) 

1 
1 
1 
1 
1 
1 
1 



uf)zfil-u<f)p^m + {l 
vf)zfil-u^)p,,iQ) + il 

.B\^C 



ut)zfil-vf)p,m 



vf)zfil-vr)p,siO), 



C'\ 



uf)(i - ur)z^'p^,io) + (1 - o(i - vf)z^psm 



vf){l 



P56(0) + (1 



vf){l 



v^)z^pjsiO). (A.IO) 
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